New Phytologist

Plant cold acclimation emerges from coordinated adjustments in photosynthesis, primary metabolism, and intracellular carbon allocation. Yet, the regulatory role of subcellular metabolite compartmentation in natural variation of cold acclimation remains insufficiently understood. Here, we investigated four Arabidopsis thaliana accessions grown either individually or in bulk to determine how growth configuration and genotype shape the metabolism of sugars and organic acids during cold exposure. Using non-aqueous fractionation, we quantified plastidial, cytosolic, and vacuolar sugar pools alongside whole-cell carbohydrates, organic acids, enzyme activities, photosynthetic parameters, and stress markers. A neural-network classifier revealed that subcellular sugar distribution together with sugar amounts and organic acids provided the strongest discriminatory power among accessions, surpassing photosynthetic traits and enzyme activities. Our findings demonstrate that natural variation in cold acclimation is strongly determined by genotype-specific subcellular metabolite architectures, and that the cultivation strategy modulates these intracellular signatures. We conclude that subcellular compartmentation of metabolites represents a cellular control layer for natural variation of cold acclimation and resilience in Arabidopsis thaliana.

Extrafloral nectaries (EFNs) are a widespread plant defense mutualism trait and are highly convergent, appearing in hundreds of plant lineages worldwide. Here we investigate a report of possible EFNs in Erythranthe angulosa, a recently-described California wildflower. We integrate field observations, insect bioassays, an induction experiment, and microscopy to test for signatures of EFN function, finding no evidence that the distinctive axillary swellings produced by E. angulosa function as EFNs. We also uncovered two distinct morphs at the type locality of E. angulosa that diverge in the number of axillary swellings produced, as well as other shoot architecture traits such as stem thickness, leaf size, and branch number. Although the axillary swellings appear to not function as EFNs, they remain a compelling morphological variant within the yellow monkeyflowers that may perform storage or another unknown function.

Accelerating photoprotective regulation to improve carbon assimilation is a promising strategy to increase crop productivity. Although rapid non-photochemical quenching (NPQ) relaxation has been validated as a target through metabolic engineering, it remains unclear whether conventional breeding has improved this trait. Here, we investigated whether more than a century of soybean breeding enhanced NPQ relaxation alongside light-saturated carbon assimilation and seed traits. We evaluated a historical panel of 24 soybean genotypes across vegetative and reproductive developmental stages by integrating NPQ relaxation, gas exchange parameters, xanthophyll-cycle pigment profiles, expression of key photoprotective genes (VDE, PsbS, and ZEP), seed number and seed weight. NPQ relaxation parameters were not consistently associated with genotype release year, seed number, or seed weight at either developmental stage. The only exception was the amplitude of the rapidly relaxing NPQ component (AqE), which was negatively correlated with all three variables during the reproductive stage. In contrast, genotype release year was positively associated with maximum net CO2 assimilation rate (Amax), maximum carboxylation rate of Rubisco (Vcmax), maximum electron transport rate (Jmax), seed number, and seed weight, while Amax and Vcmax were positively correlated with seed number and seed weight. These findings indicate that the greater photosynthetic capacity of modern genotypes was not accompanied by faster photoprotective response. Thus, photoprotective regulation has not kept pace with gains in photosynthetic capacity under field conditions. We conclude that rapid NPQ relaxation remains an important target for synchronizing photoprotection with the high photosynthetic capacity of modern soybean lines.

Anterograde and retrograde signalling establish bidirectional communication between the nucleus and chloroplasts. Retrograde signals from chloroplasts regulate nuclear gene expression while anterograde signals from the nucleus coordinate chloroplast development and maintain cellular homeostasis. How this bidirectional signalling framework extends beyond locus-specific regulation to shape the global spatial organization of nuclear chromatin across tissues remains unclear. Although anaesthesia can alter chromatin organisation, the role of chloroplast dysfunction in these changes remains unclear. Here, we investigate how chloroplast dysfunction and anaesthesia influence euchromatin and heterochromatin organisation in Solanum lycopersicum seedlings across tissues with contrasting photosynthetic competence. Using confocal and super-resolution radial fluctuation (SRRF) imaging with quantitative multiparameter analysis, we identify distinct, tissue-specific chromatin responses to chloroplast disruption and anaesthesia. Notably, anaesthesia induces distinct spatial chromatin changes across tissues that are independent of chloroplast dysfunction, suggesting a direct nuclear response to anaesthesia rather than a chloroplast-mediated retrograde effect. These findings highlight chromatin topology as a potential quantitative biomarker of cellular disruption and provide a framework for investigating anterograde chloroplast-nucleus coordination and stress-responsive nuclear organisation in plants.

Perennial grains represent a promising frontier for sustainable agriculture, but breeding progress is constrained by the accessibility of genotyping and the difficulty of evaluating complex traits expressed for multiple years after establishment across heterogeneous environments. Phenomic selection may help address these challenges by using inexpensive, scalable, high-dimensional phenotypes collected early in development, although the robustness of such predictions across breeding cycles remains uncertain. Here, we compared genomic selection and phenomic selection across two breeding cycles of Thinopyrum intermedium (intermediate wheatgrass; IWG; Kernza(R)), comprising approximately 2,280 individuals from maternal half-sib families evaluated across multiple field sites and years. We constructed relationship matrices from genomic markers and early-life stage phenomic data, including seed and leaf color (HSV), CropReporter multispectral reflectance and indices, and cycle-specific hyperspectral reflectance sensors. Genomic models provided the strongest predictions on average across all field traits in both cycles. Among phenomic predictors, leaf HSV was consistently the most informative, whereas CropReporter and hyperspectral data showed lower and more trait-dependent performance and seed HSV provided little predictive value. Genomic, leaf HSV, and CropReporter models transferred across breeding cycles with little apparent loss of predictive ability relative to within-cycle validation, demonstrating that their predictive signals were not restricted to a single breeding cycle. Early-life stage leaf HSV emerged as a practical, accessible tool for germplasm thinning and early-stage prioritization in perennial breeding programs. Despite limited similarity among relationship matrices, multi-relationship-matrix models rarely improved prediction beyond the stronger constituent single-relationship-matrix model. Together, these results show that early-life stage phenomic data provide reproducible information about agronomic performance expressed years later, but that predictor complexity and data integration do not guarantee improved prediction.

All living organisms rely on the movement of ions across cell membranes as the fundamental physical basis of their internal energy and signaling, and plants are no exception. Plants perceive, integrate, and respond to environmental stimuli through electrical signals, classified as action, variation, and system potentials, that are coupled with calcium waves, reactive oxygen species, and hydraulic and hormonal changes to coordinate whole-organism responses despite the absence of a nervous system. Yet most studies characterize these signals using a single feature, such as amplitude or spike duration, in a single tissue, an approach that cannot establish how such signals correspond to the underlying ionic activity, mobility, and structural complexity of the signaling environment, or how this correspondence varies across organs. Here, we correlate plant bioelectrical signals with potential ionic energy flow using a multi-domain framework, combining discrete spike events, continuous waveform properties, spectral composition, and signal complexity applied to leaf, stem, and root recordings from tomato (Solanum lycopersicum) exposed to different stimulus. Electrical activity with increased stimulus strength, likely reflecting increased ionic flow, with the root showing the largest response. This suggests plant electrical signaling works as a distributed, ion-based information system, useful for stress monitoring and bio-inspired sensor design.

Antibiotic-induced disruption of plant-associated microbiomes has the potential to alter host health beyond the directly exposed generation, yet whether the effects of dysbiosis are transmitted through the seed microbiome remains unknown. Here, we investigated the intergenerational impacts of streptomycin-induced dysbiosis in tomato (Solanum lycopersicum), demonstrated that seed microbiome transfer (SMT) restores progeny microbiome function and disease resistance, and characterized the underlying physiological and genetic mechanisms. Parental streptomycin exposure altered the composition of progeny rhizosphere bacterial communities, reduced expression of defense-associated genes, and increased susceptibility to Xanthomonas perforans. Suppression of immune gene expression was strongly associated with increased disease severity, indicating that parental dysbiosis impaired progeny plants ability to mount effective immune responses. Transfer of the seed microbiome from healthy plant donors partially restored rhizosphere community composition, reduced disease severity and recovered defense gene expression of three genes. Together, our findings demonstrated that antibiotic exposure microbiome disturbance generates intergenerational legacy effects that influence plant immunity and disease susceptibility and seed microbiome transfer can counteract this dysbiosis across generations.

Water availability is critical for plants and their microbial communities, including pathogens. The plant pathogen Phytophthora cinnamomi persists in soils with fluctuating moisture, yet cellular responses to water limitation remain poorly understood in Phytophthora and oomycetes more broadly. Although we recently characterized the proteomic response of P. cinnamomi to NaCl-induced osmotic and ionic stress, its response to PEG-mediated water limitation remains poorly understood, leaving a critical gap in our understanding of drought-relevant stress adaptation. Here, we quantified mycelial growth and profiled time-resolved proteome dynamics of P. cinnamomi during polyethylene glycol (PEG-3350)-treatment, simulating moderate water limiting conditions. Treatment with 5% PEG-3350 enhanced radial mycelial growth relative to controls, with no early growth inhibition observed. Label-free proteomics identified 1,097 protein groups, with 880 proteins shared between conditions and an asymmetric abundance profile dominated by decreasing protein abundance over time. Only a small subset of proteins increased, mainly enzymes involved in redox buffering (e.g., thioredoxin and glutaredoxin-like proteins) and mitochondrial/metabolic regulation (e.g., alternative oxidase) and mitochondrial/metabolic regulation. Hierarchical clustering revealed a potential three-phase temporal program: early translational and regulatory remodeling (1-6 HPT), sustained metabolic adjustment (6-12 HPT), and delayed engagement of redox and proteostasis functions (12-24 HPT). Network analysis demonstrated that redox-associated function was integrated throughout this adaptation, with individual clusters further specialized by cofactor preference (NADP- versus NAD-dependent enzymes) and distinct metabolic roles (malate dehydrogenase, CoA-ligase activity). This coordinated, multi-phase reorganization sustained mycelial growth despite moderate osmotic stress, indicating that P. cinnamomi employs active proteomic adaptation rather than passive stress tolerance. These findings reveal the cellular mechanisms underlying drought persistence in this invasive pathogen and suggest molecular targets for disease management under water-limited conditions.

MicroRNAs (miRNAs) are key post-transcriptional regulators of development, stress response, and secondary cell wall formation in woody plants, yet annotations for Eucalyptus grandis, the world's most widely planted hardwood, remain fragmented across studies using incompatible discovery pipelines and filtering criteria. Here we present the Eucalyptus MicroRNA Archive (EMA), a curated, locus-resolved database integrating three independent small RNA sequencing datasets spanning vegetative tissue, somatic embryogenesis, and mechanically induced tension wood formation. Applying annotation criteria aligned with current plant miRNA standards, EMA catalogs 99 curated miRNAs (31 previously described, 68 novel) organized into 34 family-level groupings under a three-tier confidence system, known-reference-supported, multi-study replicated, or single-study, that preserves study-of-origin and sample-level evidence for every entry. Cross-study comparison showed that only 9 of 99 entries (9.1%) were independently supported by all three datasets, supporting an evidence-tiered rather than binary annotation scheme. Target prediction against the E. grandis transcriptome yielded 1,773 miRNA-target interactions spanning 764 loci, integrated into a combined miRNA-target and protein-protein interaction network. This network resolved into functionally coherent, mutually isolated clusters, including an miR482-associated NBS-LRR/TIR disease-resistance hub with a substantial translational-repression component, alongside modules enriched for ribosome biogenesis and translation, DNA replication, and nitrogen and carbohydrate metabolism. EMA is publicly accessible through an interactive web dashboard, with all curated data, source code, and analysis scripts openly available, providing a reproducible, extensible framework for E. grandis miRNA research and a template for similarly structured resources in other non-model woody species.

Although the biosynthetic pathways of benzylisoquinoline alkaloids (BIAs) have been extensively investigated in several plant species, their transcriptional regulatory mechanisms remain only partially understood. Jasmonate (JA)-responsive group IX APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors (TFs) are well-known regulators of specialized plant metabolism, including the biosynthesis of various alkaloids. However, their specific roles in BIA biosynthesis remain largely elusive. Here, we isolated five novel group IX AP2/ERF TFs, designated Benzylisoquinoline alkaloid Jasmonate-responsive AP2/ERF (BJE1-5), from Coptis japonica. Phylogenetic analysis revealed that Benzylisoquinoline alkaloid Jasmonate-responsive AP2/ERF (BJE) proteins belong to subclades distinct from group IXa, which contains well-known AP2/ERF TFs involved in alkaloid biosynthesis. Transient expression analyses in C. japonica protoplasts demonstrated that certain BJEs, particularly CjBJE3 and CjBJE5, positively regulated BIA biosynthetic genes through a mutual regulatory network among BJE members. Moreover, CjBJE3 expression was regulated by CjbHLH1, a unique-type basic helix-loop-helix (bHLH) TF specific to BIA-producing plants. Furthermore, heterologous expression of CjBJE3 and CjBJE5 in cultured Eschscholzia californica cells significantly enhanced the overall BIA production, particularly by increasing end-product benzophenanthridine BIAs, highlighting several uncharacterized biosynthetic genes clustered in the genome. Our findings suggest that BIA-producing species have developed a specific regulatory network comprised of CjbHLH1 and BJE TFs, providing valuable clues for identifying novel biosynthetic enzymes.

Aegilops peregrina is a wild allotetraploid wheat wild relative and an important source of genetic diversity for stress tolerance and agronomic traits. Here, we report a subgenome-resolved, chromosome-scale reference genome assembly of a drought tolerant and stem rust resistant Ae. peregrina accession PI 604178 generated using PacBio HiFi and Hi-C sequencing. The 10.13 Gb assembly contains 98.81% of sequence anchored to 14 pseudomolecules representing the seven S and seven U chromosomes, with contig and scaffold N50 values of 25.84 and 746.48 Mb, respectively. The assembly achieved a consensus quality value of 74.61, 97.83% k-mers completeness, and 99.9% BUSCO completeness. LTR Assembly Index values of 20.43 and 18.79 for the S and U subgenomes, respectively, further supported high continuity across repeat-rich regions. Repetitive elements comprise 85.93% of chromosome-anchored assembly. We annotated 59,910 high-confidence protein-coding genes, with comparable gene representation across the two subgenomes. This reference genome provides a high-quality genomic framework for comparative analyses, characterization of important loci regulating agronomic and resilience related traits, and sequence-guided exploitation of Ae. peregrina allelic diversity for wheat improvement.

Blumeria graminis f. sp. avenae (Bga), the causal agent of oat powdery mildew, is one of the most host-specialized members of the B. graminis species complex. Despite its agricultural importance, the lack of a high-quality reference genome has limited studies of host specialization, virulence evolution and comparative genomics in this pathogen. Here, we generated the first chromosome-scale genome assembly of Bga using an integrative approach combining long- and short-read sequencing, Hi-C scaffolding and transcriptome data. The Bga genome exhibits hallmark features of powdery mildew fungi, including extensive repeat content and low gene density. Comparative analyses revealed that genome expansion is primarily associated with historical transposable element proliferation rather than recent transpositional activity. Genome organization is consistent with a functionally stratified "one-speed" model, in which genes associated with pathogenicity, including predicted effectors and infection-responsive genes, are preferentially located in transposable element-rich regions characterized by reduced synteny conservation and extended intergenic spaces. In contrast, conserved genes are concentrated in compact genomic regions and maintain strong syntenic conservation across cereal-infecting formae speciales. Hi-C analyses demonstrated a highly structured chromatin architecture and revealed genome organization patterns associated with infection-related gene expression. Comparative genomic analyses indicated that host specialization in Bga is driven by localized diversification of a relatively small subset of genes rather than large-scale genome restructuring. These results provide the first high-quality genomic resource for Bga and offer new insights into the evolutionary mechanisms underlying host specialization in powdery mildew fungi.

Zinc is an essential micronutrient at low concentrations, yet it becomes toxic at slightly higher ones. This is exploited by plants as an effective defensive strategy. However, the molecular components that are involved zinc-mediated immunity remain poorly defined. Here, we show that mixed-linked {beta}-1,3/1,4-glucans naturally occurring in microbial and grass cell walls and used as an agrobiological solution, trigger zinc accumulation in the Arabidopsis apoplast and upregulate the expression of the zinc transporters HMA2 and HMA4. This response occurs independently of salicylic acid, jasmonic acid and ethylene-mediated signalling pathways, but it requires the LysM receptor kinases CERK1, LYK4 and LYK5, indicating a specific pattern triggered immunity-associated mechanism. We further demonstrate that hma2hma4 mutants display constitutive activation of a broad set of defence-related genes, yet this transcriptional reprogramming is insufficient to confer resistance against the necrotrophic fungus Plectosphaerella cucumerina BMM. Moreover, metabolomic profiling highlights the contribution of specialized metabolites to this defective defence output. Altogether, our findings reveal that zinc-mediated toxicity constitutes a defence mechanism integrated into the immune response triggered by specific microbial or damage associated molecular patterns.

1. Forest canopies harbour strong microclimatic gradients that shape plant performance, species interactions and ecosystem processes. Yet, despite renewed interest sparked by global change, forest canopies remain difficult-to-access experimental spaces. 2. With the goal to expand access to tree canopies as experimental arenas, we designed, built, and tested TreeTOP, a standardized experimental platform that opens canopy space for manipulative ecological experiments, specifically with potted plants. TreeTOP features lightweight aluminum frames placed in mature tree canopies non-invasively, allowing potted plants to be placed in three different heights, ground level, shade canopy, and sun canopy. 3. We implemented TreeTOP using two contrasting infrastructure concepts to demonstrate its applicability in both highly equipped canopy research facilities and forests without permanent canopy infrastructure. One installation relied on a canopy crane, grid power and fully automated irrigation, whereas the second was built by certified tree climbers and was equipped with an autonomous solar-powered, battery-operated irrigation system. At both sites, environmental sensor networks monitor the experiment. 4. TreeTOP successfully reproduced characteristic canopy microclimatic gradients, including increasing light availability, daytime air temperatures and thermal extremes with canopy height. Despite differing infrastructures, both implementations generated comparable microclimatic patterns, demonstrating that standardized canopy experiments are feasible in forests with or without permanent canopy access. By opening canopy space for manipulative experiments, TreeTOP provides a transferable framework for investigating plant performance, phenology, species interactions and microbiome assembly under realistic forest conditions.

Brettanomyces bruxellensis is a yeast species associated with diverse fermentation environments and characterized by extensive genetic diversity, including diploid, autotriploid, and allotriploid lineages resulting from independent hybridization events. These lineages are associated with distinct ecological niches and provide a framework for studying metabolic trait evolution in complex genomes. Nitrate assimilation is a relatively uncommon trait among yeasts and has been reported in B. bruxellensis, but its distribution and evolutionary history within the species remain poorly understood. Here, we combined phenotypic characterization of 151 strains with genomic analyses of 946 whole-genome sequences to investigate nitrate assimilation. Growth assays revealed that nitrate assimilation is widespread but unevenly distributed across genetic lineages, with some populations largely retaining the trait whereas others have frequently lost it. Genomic analyses identified extensive variation affecting the nitrate assimilation gene cluster composed of YNR1, YNI1, and YNT1. Nitrate assimilation was strongly associated with both gene copy number and predicted gene functionality, with nitrate-assimilating strains generally carrying more functional copies of the cluster. Leveraging the complex genomic architecture of the species, we independently analyzed primary and acquired genomes in allotriploid lineages and uncovered contrasting evolutionary trajectories following hybridization. While nitrate assimilation genes were generally maintained in primary genomes, acquired genomes showed a higher prevalence of gene loss and predicted loss-of-function variants, revealing asymmetric dynamics between subgenomes. Altogether, our results suggest that nitrate assimilation represents an ancestral trait that has been differentially maintained across B. bruxellensis lineages through a combination of copy number variation, gene degeneration, and genome-specific evolutionary dynamics. These findings provide new insights into how genome architecture and polyploid evolution shape the maintenance and loss of metabolic traits in an industrially relevant yeast species.

Water-use efficiency (WUE), the ratio of accumulated plant biomass to water lost through transpiration has conventionally been determined using a destructive single-point measurement. Recent advances in high-throughput phenotyping now enable repeated, non-destructive estimation of biomass and WUE. However, these digital measurements must be statistically validated against conventional destructive methods to validate their use as reliable proxies. Therefore, we compared digital biomass determined point clouds produced from multispectral camera scanners with destructive harvests across eight harvests using Samsun tobacco grown under both drought and high-water conditions. WUE efficiency, calculated using the digital biomass estimated from a point cloud and gravimetric water use determinations, were compared to destructive harvest determinations. The coefficient of variation (CV) showed there were no significant differences in digital and destructive measurements for either biomass or WUE. Indicating that digital measurements can be used in place of destructive measurements. Drought plants used significantly less water and were significantly smaller than high-water plants from Harvests 4 through 8. However, there were no significant differences in the ratio of evapotranspiration to leaf area or WUE, indicating that drought plants were simply smaller and used less water than the high-water plants. This work validates that estimating plant biomass from a digital point coupled with continuous gravimetric determination of water use provides a reliable nondestructive measure of WUE in high-throughput measurements across the full plant life cycle.

Drug resistance emerges repeatedly in outbreaks of Candida fungal pathogens, but little is known about its origins or persistence. Here, we investigated the evolutionary processes shaping echinocandin resistance in Candida auris, a globally emerging and predominantly clonal fungal pathogen. Genome-wide association across over 600 isolates identified mutations in the {beta}-1,3-glucan synthase gene FKS1 as the most significant driver of resistance to an echinocandin drug. Ancestral reconstruction of this population traced shared resistance mutations among small groups typically consisting of 2-3 closely related isolates, but clusters could include up to 16 isolates. Nearly all resistant clusters consisted of isolates collected in the same year and region, consistent with local transmission. To further examine population-level selection, we measured adaptive signatures in FKS1 and the highly diverged paralog FKS2 across 22,000 genomes. This revealed excess nonsynonymous polymorphisms in FKS1, primarily due to independent, recurrent mutations at resistance hotspots, consistent with parallel evolution and incomplete fixation of adaptive alleles. In FKS2, there is no evidence of hotspots and little support for diversifying selection. Together, these results indicate that resistance mutations emerge under strong genetic constraint, with adaptation restricted to only one FKS homolog and predominantly at mutational hotspots.

The integration of digital technologies for high-throughput field phenotyping is critical for accelerating crop improvement in agriculture. However, extracting traits from remote sensing data remains constrained by fragmented workflows, manual intervention, and limited interoperability among existing tools, resulting in delays that hinder timely biological insight and decision-making. To address these challenges, we present PhenoStream (Phenotyping Streaming), a scalable, end-to-end cyberinfrastructure designed to automate the full lifecycle of aerial imagery-based phenotyping, from data acquisition to plot- and genotype-level inference. The framework integrates automated data ingestion from distributed field sites, geospatial processing, and AI-enabled trait extraction within a unified, user-accessible graphical interface. Its modular and extensible architecture supports adaptable trait modeling and seamless integration of new data sources, enabling deployment across diverse crops, environments, and experimental designs. We demonstrate the system across a large multi-location field trial network of bioenergy crops, where it enables high-throughput characterization of spatiotemporal growth dynamics, genotype-by-environment (GxE) interactions, and predictive modeling of key agronomic traits. By significantly reducing processing latency and manual effort, the platform facilitates near-real-time analysis and reproducible workflows. This work establishes a generalizable and scalable pathway for operationalizing very-high-spatial resolution aerial phenotyping in agricultural research. By bridging data acquisition and analytics, the end-to-end cyberinfrastructure provides a foundation for integrating heterogeneous and unstructured data streams--including remote sensing, environmental, and management data--toward data-driven decision making in agriculture.

High-phenolic plant genera present a major technical limitation in genomic research. Standard extraction approaches that perform reliably across diverse flora often perform poorly when applied to recalcitrant taxa, producing low DNA yield and integrity incompatible with sequencing requirements. The genus Anigozanthos (Kangaroo paws) from the family Haemodoraceae exemplifies this problem. We identified key physicochemical factors governing extraction failure in this genus and resolved them through targeted modifications to lysis chemistry and contaminant management. The resulting protocol achieved a near threefold improvement in DNA purity, substantially reducing contaminant carry over and consistently yielded high-integrity, long DNA fragments (DIN > 7) across a diverse sample set spanning cultivated and wild material across four diverse genera of Haemodoraceae. We also tested a straightforward purity assessment framework that can be implemented in any standard molecular laboratory, enabling rapid pre-submission quality assessment without the need for specialised equipment. Together these advances open a practical path to genomic characterisation of Anigozanthos that establishes a transferable model for genomic research across Australia ' s chemically complex native flora.

Oil palm (Elaeis) is the primary source of global vegetable oil. Interspecific hybrids of Elaeis exhibit pronounced heterosis by integrating two distinct subgenomes into a single nucleus, effectively combining the high yield of African oil palm (E. guineensis) with the high unsaturated fatty acid content and disease resistance of American oil palm (E. oleifera). However, the genetic basis underlying heterosis is still unclear. Here, we combine phased genome assembly, comparative genomics, evolutionary genomics and haplotype-aware transcriptomics to unravel the genetic architecture of heterosis of hybrid oil palm. We assemble the highly heterozygous F1 genome ('Reyou 40', 3.75% heterozygosity) into a complete 1.73 Gb T2T haplotype (HapG) and a 1.84 Gb near-T2T haplotype (HapO with17 gaps). Despite 91.56% sequence identity, HapG and HapO diverged in LTR-RT occurrence and PAV affected genes, showing complementary biases in lipid metabolism and stress responses, respectively. Evolutionary genomics revealed that ancient WGDs preserved the palm family. Whereas lineage-specific lipid-related gene expansions in oil palm. Six ancient introgressed regions (~64 Mb) in HapG were reshaped by transposable elements and tandem duplication, showing an enrichment of genes related to resistance and lipid metabolism. Transcriptomically, 82.2% of allelic gene pairs maintained balanced expression, accompanied by parental functional complementarity and dosage buffering, revealing a potential regulatory basis for coordinating parental genetic differences in the hybrid genome. These haplotype-resolved genomic resources offer vital targets for understanding heterosis and accelerating oil palm molecular breeding.